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Chapman's model for ozone concentration is studied. In this nonlinear 
model, the photodissociation coefficients for O2 and O3 are time-depending 
due to earth-rotation. From the Kapitsa's method, valid in the high fre- 
quency limit, we find the criterion for the existence of equilibrium solutions. 
These solutions are depending on the frequency, and require a rotation period 
T which satisfies T < Ti or T > T2. Where the critical periods Ti and T2, 
with T2 > Ti, are a function of the parameters of the system (reaction rates 
and photodissociation coefficients). Conjectures respect to the retardation of 
the earth's rotation, due to friction, suggest that the criterion was not even 
verified in the atmospheric past. 
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1.- Introduction 



The dynamics of the ozone layer in the atmosphere has different basic 
process hke: chemical reactions, photochemical reactions and transport (dif- 
fusion, convection, etc.). In a general point of view, this dynamics is complex 
and requires some approximations to be studied. In this sense, we consider a 
photochemical model proposed by S. Chapman (Brasseur, 1986; Chapman, 
1930; Wayne, 1991). This model considers a set of reactions between the 
oxygen components. Explicitly, 



In the reaction of the ozone production Rl, M denotes any atmospheric 
element acting as catalyzer. The reaction R2 denotes the loss of oxygen O 
and ozone O3 producing molecular oxygen 02- R3 and R4 correspond to 
photochemical destruction process related to the solar radiation (symbolized 
by hu). 



The time evolution equations for the constituents consider the above re- 
actions for variable concentration assuming the concentration of O2 being 
stationary. Let X and Y be the concentration of O and O3 respectively. 
Then, the Chapman's model (see for instance Brasseur, 1986; Montecinos, 
1998; 1999) considers the time evolution equations for the concentrations 
given by 
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R2) 
R3) 
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O + O2 + M ^03 + M, 
O + O^^ 2O2, 

02 + hu^O + O, 

03 + hiy^0 + O2. 
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Ji + J2Y - {h + k2Y) X, 



(1) 
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dY_ 
It 



hX - ( J2 + k2X) Y, 



(2) 



where, on the right hand, the positive terms are production rates, and the 
negative ones are loss rates. In the nonhnear system (1,2), the quantities 
Ji and J2 are related to the reactions R3 and R4 and correspond to the 
photodissociation of O2 and O3 respectively. The important fact is that they 
dependent on the sun's radiation, and then, they are periodic in time with a 
period T — 2A hours. In this paper, and by simplicity, we assume 



where uj = 2tt/T , and J° are positive constants. On the other hand, the 
positive constants ki and k2 in (1,2), are temperature dependent. Also they 
are dependent on the O2 concentration, and related to the reaction velocity 
in Rl and R2 respectively (DeMore, 1994). 



In this paper we propose an analytical study of the nonlinear model (1- 
3). In a general point of view, the study of this system is a difficult task, 
nevertheless, some interesting results can be find in the high frequency (a;) 
limit. In fact, we use a method proposed originally by Kapitsa (Landau, 
1982) for a mechanical particle in a field with rapid temporal oscillations in 
the parameters and nonlinear terms. We find explicitly the solutions of the 
systems (1-3) in the high frequency regime (19,20). They are the equilibrium 
solution for the time-averaged concentration. Calculations tell us that no 
solution exist in some frequency range. The study of the behavior of this 
model, for different frequencies, has a physical interest because earth period 
of rotation changes with the geological age. It is a known fact that the 
rotation of earth is gradually slowing down by friction. In fact, the estimated 
rotation velocity diminishes by 4.4 hours every billon of years (Shu, 1982). 



For explicit calculations, we assume the following values for the parame- 
ters: 



Ji" ~ 10^[l/s]; ~ 10-^[l/s]; fci ~ 10[l/s]; ^2 ~ 2.5 x 10~^^[l/s], (4) 
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J.(i) = J°(l -cosa;i), i = l,2 



(3) 



corresponding to the values for the ozone layer altitude at, more or less, 35 
km. 



Finally we note that the case of small frequency (a; — > 0) can be solved. 
In fact, here the parameters Ji{t) evolve slowly with time and then, they can 
be assumed as constant in the integration process of (1,2). So the solutions 
corresponding to the 'fixed point' ~ ^ ~ 0, for definitions see (Seydel, 
1988; Wio, 1997)) are 



4JJ \| 2k2Jf 

Namely, in this approximation, the variable Y is constant and X varies 
linearly with the dissociation coefficient. Moreover, the solutions (5,6) are 
consistent with the numerical solution in the stratosphere (Fabian, 1982; 
Montecinos, 1996). We note that, from (5), it is easy to show that the 
variable Y satisfies F < In fact, it is a general bound for the solution of 
the systems (1,2) (Montecinos, 2000). 



2.- The method of Kapitsa 

As said before, we shall study the system (1-3) in the high frequency 
regime. High frequency means here small period of oscillations T respect to 
the relaxation-time Tr for the slow variables. 



Assume the separation of the concentrations X and y in a slow temporal 
variation {x and y) and other fast {e and r), respectively). Namely, 

X^x + s; Y^y + T], (7) 
where the fast variables are periodic with temporal average zero, namely. 



4 



(^)t = iv)T = 0. (8) 
The equation (1) can be re-written like 

dx d,£ 

+ ^ = + - ih + k2y) X - ( J° + J^y) cosu;t+ 

{J° — J2 cosLut — k2x) T] — {ki + k2y) £ — k2S'r), (9) 
and equation (2) becomes 

^ + ^ = kix — ( J2 + k2x) y + J2y cosu)t— 
dt dt 

(J| + k2X — J2 cosa;t) rj + {ki — k2y) e — k2eri. (10) 

On the other hand, the fast variables are only related to rapid oscillation 
(Landau, 1982). In this way, from the above expression (9,10), they are 
assumed to be a solution to the differential equations: 

ds dn 

— = - ( J° + J-y) COS cut; = J°y cos ut. (11) 

At this point a remark becomes necessary. The expression (7), comple- 
mented with the above equations (11), defines a change of variables without 
approximations. Nevertheless, the differential equations (11) are suggested 
by the direct oscillatory term in (9) and (10). Kapitsa's method consider the 
equations for ^ and ^ as approximated. 



In one period, the slow variables are essentially constants and the fast 
have zero average (8), then, the time average of the equation (9) becomes. 

dx 

— = Ji + J2I/ - (^"1 + ^2y) X - J2 {f] cosa;t)y - /c2 {er])T > (12) 
and, for equation (10), we obtain 

-J- = kix — ( J2 + k2x) y + J2 {v cosujt)j, — k2 {£fl)j, . (13) 
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Nevertheless, since equations (11) can be solved exactly, 



s{t) 



1 



JO 

(J° + J2?/) sinwi; r]{t) = —ysmcijt, 



(14) 



the evolution equations for the slow variables become 



dx 



J1 + J°y - {h + k2y) X + 



(15) 



and 



dt 



kix — ( J2 + k2x) y + 



J!ik2 



y{Jt + J2y)- 



(16) 



This set of equations are the basis for our analytical results. They are 
restricted to the high frequency approximation. This approximation becomes 
given by the 'expansion' in l/o;^ related to the last term in (15,16). Remark 
that it is an autonomous nonlinear systems and then, without the explicit 
temporal dependence. This transformation, from a set of equations with 
time-periodic parameters, to other autonomous, is related to the Kapitsa 
original ideas (Landau, 1982). 



In a general frame of work, the system (15,16) is complex. Moreover, 
the approximation of high frequency is valid when the relaxation time Tr, 
of the equations (15,16), is bigger than 27r/cj. This comparison is a difficult 
task, nevertheless, the case J| = can be solved exactly to estimate the 
validity of the approximation. It corresponds formally to eliminate the disso- 
ciation of 0,3. The asymptotic solution of the non- autonomous system (1-3) 
is (Montecinos, 2000). 



where the phase is given by the relation: tan0 = uj/2ki. On the other 
hand, combining the equations (15,16) with (7), in the high frequency ap- 
proximation we found that the equilibrium solution given by the Kapitsa's 
method is: 




(17) 





(18) 
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It is direct to show that the exact solution (17) reduces to (18) in the high 
frequency limit. Moreover, the relaxation time Tr can be calculated here 
analytically. It is given by Tr = l/2A;i. So, we expect that the high frequency 
approximation (15,16) is valid when Airki <^ cu. 



3.- Existence of equilibrium solutions 

In this section we are concerned with the fixed point solution (Seydel, 
1988; Wio, 1997) of the autonomous set (15,16). This system have an equi- 
librium point (xo,yo), defined by ^ = ^ = 0, and given by the solution of 
the equations: 

k,J^ (2-^]yl + hr, (l -^]yo- hr, = 0, (19) 



and 



+ .7?y„) (l + ^yo) 



Equations (19) and (20) define the homogeneous equilibrium solution of 
the autonomous system (15) and (16) and then, with (7) and (14), we have 
the solution of the systems (1,2) in the high frequency regime. 



The existence of real solutions, for the second degree equation (19), re- 
quires the inequality 

1 4 / k,Jt\ 1 8 / fe^^i°Uo (21) 

which corresponds to an inequality of second degree for 1/a;^. Namely, there 
is no equilibrium solution of (15,16) if and only if. 



2 1 2 k2J? 12 12 k2J? 
\ \ 1 — < — < \ \ \ 1 — 
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If we assume the condition 



k2J\ 




(23) 



valid for the parameters (4) of section 1, the inequahty (22) can be re- written 
for the period T. In fact, there is no equilibrium solution of the system (15,16) 
when 



4.- Earth's slowing rotation and the existence 
of solution 

It is interesting that the inequality (24) gives a region were no solution 
exist. Here we must take care because no oscillating solution like (7) exist. 
In fact, (24) splits the a;— space parameter in three regions: (i) The region 
defined by T < Ti where a real positive solution (|/o > 0) of equation (15) 
exist, with a negative one < 0) . (ii) The region defined by (24), where 
no solution of (15) exist, (iii) The region defined by T > T2 where solutions 
are negative {ijo < 0). Prom (5,6), we known that in the slow frequency limit 
(region (iii)) a real positive solution exist. Then. Kapitsa's method does 
not work well in this region, nevertheless, at least it says that an oscillating 
solution exist. 



At this point we can formulate the following question: since the earth- 
rotation has diminished by friction, how has the change in rotation affected 



Ti < T < T2 {no - solution), 



(24) 



where 




(25) 
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the existence of the ozone layer ?. This question seems appropriate because 
the Kapitsa's method tells us that the frequency of rotation and the ozone 
concentration are related. Using the parameter values (4), of section 1, we 
can estimate the critical period (24) : Ti ~ 0.02 hours and T2 ~ 22 hours. Is 
interesting that the actual period T = 24 hours, is in the region of permitted 
solution (T > T2). Moreover this is suggestive: from the retardation of 
earth rotation velocity data (4.4 hours/billon of years, (Shu, 1982)), a simple 
calculation tells us that before ~ 0.46 billons years no solution existed 
because we were in the region (24). This is a surprising estimation if we 
consider that actually the ozone layer is believed to have been in existence 
0.7 billon years (Graedel, 1993). 



5.- Conclusions 

We have considered the Chapman's model for ozone production (1,2). In 
this nonlinear model, the parameters related to photodissociation are periodic 
in time (3). We were interested at the analytical study of this model by 
using the high frequency approximation, due to Kapitsa. Namely, we have 
considered the autonomous system (15,16), depending on the frequency, for 
the averaged variable concentrations. The existence of equilibrium solutions 
(fixed points (19,20)) is depending on the frequency. In fact, there are two 
critical period Ti and T2 so that for Ti < T < T2 there is no equihbrium 
solution (24). 



The values for the parameters (4) , in section 1 , give the condition of no- 
existence (24): 0.02 hours < T < 22 hours, and then compatible with the 
actual earth's period of rotation, and existence of the ozone layer. Moreover, 
considering the earth's slowing rotation motion due to friction (4.4 hours 
every billon of years, (Shu, 1982)), we estimate that the ozone existence 
condition is verified after ~ 0.46 billon years (section 4). This is a good 
estimation if we consider the simplicity of the autonomous model given by 
equations (15,16). The age of the ozone layer is 0.7 billons of years (Graedel, 
1993)). 
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Before to ending a remark, equations (15) and (16) are very adequate to 
the study of diffusion process, wfiicfi was neglected in tfie original equations 
(1) and (2). In fact, because they are not time depending, when we add 
spatial diffusion terms D-^X and D-^Y, they become similar to reaction- 
diffusion-equations . 
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